Dynamical Casimir effect in a Josephson metamaterial 



Pasi Lahteenmaki 1 , G. S. Paraoanu 1 , Juha Hassel 2 , and Pertti J. Hakonen 1 
1 Low Temperature Laboratory, Aalto University, PO Box 15100, FI-00076 AALTO, Finland and 
2 VTT Technical Research Centre of Finland, PO BOX 1000, FI-020U VTT, Finland 



Vacuum modes confined into an electromagnetic 
cavity give rise to an attractive interaction between 
the opposite walls [Y[. When the distance between 
the walls is changed non-adiabatically, virtual vacuum 
modes are turned into real particles, i.e. photons are 
generated out of the vacuum [2, 3j. These effects are 
known as the static and dynamical Casimir effect, 
respectively. Here we demonstrate the dynamical 
Casimir effect using a Josephson metamaterial em- 
bedded in a microwave cavity at 5.4 GHz. We achieve 
the non-adiabatic change in the effective length of the 
cavity by flux-modulation of the SQUID-based meta- 
material, which results in a few percent variation in 
the velocity of light. We show that energy-correlated 
photons are generated from the ground state of the 
cavity and that their power spectra display a bimodal 
frequency distribution. These results are in excellent 
agreement with theoretical predictions, all the way to 
the regime where classical parametric effects cannot 
be of consequence. 

A fundamental theoretical result of modern quantum field 
theory is that the quantum vacuum is unstable under cer- 
tain external perturbations that produce otherwise no conse- 
quences in a classical treatment ,4. As a result of this insta- 
bility, virtual fluctuations populating the quantum vacuum 
are converted into real particles by the energy provided by 
the perturbation. For example, the application of an intense 
electric fields extracts electron-positron pairs from vacuum 
(Schwinger effect), the bending of space-time at in the intense 
gravitational field at event horizons results in black hole evap- 
oration (Hawking radiation), the acceleration of an observer 
in the Minkowski vacuum results in the detection of particles 
(Unruh effect), and sudden changes in the boundary condi- 
tions of electromagnetic field modes or in the equations of mo- 
tion creates photons (dynamical Casimir effect) [5]. None of 
these dynamical effects have been yet experimentally verified, 
although preliminary evidence for the analog of Hawking ra- 
diation and dynamical Casimir effect have very recently been 
reported [HE]. 

In general, the dynamical Casimir effect (DCE) is the phe- 
nomenon of particle creation from vacuum due to modulation 
of the background in which a quantum field propagates [8]. 
There are several proposals how to observe these phenomena 
using single vibrating mirrors or vibrating cavities [9l412|. In 
addition, several alternative proposals for experiments that 
could lead to a measurable flux of Casimir photons exist, for 
example sonoluminescence [13], changing the length of three- 
dimensional cavity resonators by using a laser-modulated re- 
flectivity of semiconducting layers [144416] . surface acoustic 
waves [T7], and superconducting qubits [18]. The experimen- 
tal verification of the dynamical Casimir effect in mechani- 
cally modulated systems has so far been hampered by the fast 
speeds required for the change in order to obtain a measur- 
able amount of energy [5]. A moving mirror creates photons 
out of vacuum fluctuations in pairs whose frequencies sum up 
to the pump frequency of the mirror. Typically, the photon 



production is non-negligible only when the mirror velocity 
approaches the speed of light. For example, a mirror on a 
1 GHz nanomechanical oscillator moving at an amplitude of 
1 nm would create only about 10 -9 Casimir photons/s. Re- 
cently, it was suggested that these speeds can be obtained in 
a circuit QED (quantum electrodynamics) setup in which the 
boundary condition is realized by a SQUID (superconducting 
quantum interference device) that terminates a superconduct- 
ing waveguide 20HI22]. 

In our work, we use a metamaterial consisting of an array 
of 250 SQUIDs, embedded into a 7 mm long superconducting 
coplanar waveguide. In order to enhance the photon pro- 
duction rate even further, we operate the SQUID array as a 
low Qrcs resonator (Q rcs ~ 100) by coupling the metamaterial 
through a capacitor to the external transmission line (see Figs. 
[TJi and [lb). The flux-tunable SQUIDs allow us to modulate 
the speed of light (v w 0.5co, with cq the speed in vacuum) 
and thus change the effective electrical length of the cavity 
at drive frequencies uJd on the order of 10 GHz over a rela- 
tively large range, up to 10 mm. Therefore, the field inside 
the cavity experiences an effective modulation of boundary 
conditions (see Fig. la). Moreover, the input-output theory 
[231 124] predicts that the field reflected by the cavity acquires 
a phase that depends on the detuning with respect to the 
resonant frequency of the cavity. From the point of view of 
an observer detecting the output field, this change of phase 
could have come as well from an effective mirror reflecting 
the input field or from an other appropriately chosen bound- 
ary condition. The modulation of the resonance frequency is 
thus analogous to a modulation of such a boundary condition. 
The detailed arguments of this equivalence are outlined in the 
Supplementary material. The resonance conditions facilitate 
the creation of a large amount of radiation even at temper- 
atures where the field modes are essentially unoccupied, and 
we detect a clear correlation at two frequencies symmetric 
around half the drive frequency uJd/2. We observe two-mode 
squeezed (TMS) photons and a characteristic double peak [21] 
in the squeezed output noise spectrum when the resonator is 
detuned from tdd/2. 

The samples were fabricated using a NbAlO^Nb trilayer 
process. The metamaterial cavity was cooled down to T = 50 
mK in a pulse-tube-based dilution refrigerator; T corresponds 
to a thermal occupation number of n — 0.0056 at 5.4 GHz. 
All the lines used to flux bias the sample and to feed the 
microwave signals are strongly attenuated: the microwave at- 
tenuation 7 > 160 dB in the DC-feed and, for the microwave 
cables, 7 ~ 70 dB to guarantee proper thermalization of the 
background radiation. The connection to the low noise ampli- 
fier (LNA) is isolated by two circulators in combination with 
high- and low-pass filters at the base temperature, providing 
isolation on the order of 40 dB from the 4 K stage. 

We first established that the device worked indeed as a 
cavity whose effective length is controlled by an external mag- 
netic field. Fig. [TJ; shows the variation of u) res /2Tr with mag- 
netic flux bias $. The inset displays the measured phase of the 
reflection coefficient Su which changes steeply while sweeping 
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FIG. 1: (a) Equivalent electrical and mechanical circuits: the 
modulation of the Josephson inductance in the metamaterial 
by magnetic flux $ varies X/d, the wave length with respect 
to the cavity length, which is equivalent to modulating the 
effective length d of the cavity by mechanical means. The 
coupling capacitor is equivalent to a semitransparent mirror, 
(b) Schematics of the measurement setup. The metamaterial 
sample is a 7-mm-long coplanar waveguide with 250 embed- 
ded SQUIDs, each junction having a critical current of ~ 10 
fj,A. Modulation of <E> through the SQUIDs is applied using a 
lithographically fabricated spiral coil underneath the metama- 
terial. c) Resonant frequency ui ICS /2n vs. reduced magnetic 
flux & C xt/&o without the pump signal; the DC operating point 
for DCE experiments is denoted by a green circle. The inset 
displays the measured phase of the scattering parameter Sii 
while sweeping $. The steepness of the variation in the phase 
arg(Sn) governs the effective "movement of the mirrors". 



$. The variation of arg(Sii) is a direct measure of how much 
the wavelength A ~ 4d changes with the external flux; here 
d denotes the physical length of the resonator. The change 
in X/d due to modulation is the equivalent to the " movement 
of the mirror". As shown in the Supplementary material, 
the resonance quality factor enhances the length modulation 

which becomes 5Lif = Q rcB -^- f, where I denotes 

the inductance per unit length and 51 its change due to the 
modulation. 

We also checked that our metamaterial works as a phase- 
sensitive parametric amplifier 25 - 27 in the three wave mode 
when pumped at w<j = a pump power Pd = —75 

dBm (after the low-T attenuators), the gain in the degener- 
ate mode reached 23 dB and noise temperature at 5.5 GHz 
became T n = 0.4 ± 0.2 K as deduced from the increase in 
the signal/noise ratio. Taking into account the losses in the 
circuitry of the S/N ratio measurement, this result indicates 
that the metamaterial behaves as designed. 

The existence of the Casimir photon generation was inves- 
tigated at temperatures T — 50 — 230 mK by measuring the 
noise power spectra (see Fig. 2a) when the flux was mod- 
ulated around a fixed DC bias point at twice the resonant 
frequency lj vcs /2it. To understand these results, we present 
first a simplified theoretical description of the Casimir effect 
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FIG. 2: (a) Color map of the output noise power (color scale 
in dB units with fcs x 5 K as the reference level) measured 
across the resonant frequency lj ics /2tt — 5.4 GHz while in- 
creasing the pump power Pd at oJd = 10.8 GHz. At Pd > —75 
dBm, the device becomes unstable against self-sustained os- 
cillations, (b) Cross-correlations for in-phase (I) and quadra- 
ture (Q) amplitudes (/+(*)!_(*')) and (Q+(t)Q-(t')) between 
sidebands at uj± = Ud/2±u vs. time delay t—t' with the pump 
on. The pump is at 10.8 GHz while the IQ-demodulators 
are at 5.398 GHz and 5.402 GHz. The used bandwidth is 
200 kHz and 10 6 samples were collected, (c) Covariance ma- 
trix a(A,B), measured at Pd = —74 dBm. The diagonal 
elements amount to 0.5 ■ DCE + | quanta (marked in red 
and green, respectively), while the non-zero off-diagonal ele- 
ments describe the squeezing correlations of the system, (d) 
Covariances a(x+,x~) (A), a(x+,x+) (+), and their ratio 
cr(x+,X-)/cr(x+,x+) (o) vs. modulation power Pd in the flux 
coil; the ratio and its fit are on linear scale (left), whereas the 
rest is on logarithmic scale (dB, right). The dashed and dot- 
ted curves represent theory predictions, obtained using Eqs. 
(S50-S53), (S56-S58), respectively, and the solid curve is their 
ratio. On the right scale, \ of a quantum corresponds to -6 
dB (denoted by arrow). The parameter values used in the 
evaluation are listed in the caption of Fig. 3. 



in our system. 

One can intuitively picture the dynamical Casimir effect 
as follows: mathematically, the metamaterial is equivalent 
to a gravitational pendulum with length proportional to the 
total Josephson energy Ej. The mass of the pendulum is 
the analog of the electrical capacitance of the junction and 
the angle with respect to the vertical is the analog of the 
superconducting phase difference. By modulating the effec- 
tive length at twice the natural frequency of the pendulum, 
any initial oscillation or fluctuations of the pendulum become 
parametrically amplified. However, in the absence of any ini- 
tial oscillation, changes in the length of the pendulum do not 
produce any oscillation, according to classical physics. Yet the 
quantum-mechanical prediction is different: due to the uncer- 
tainty principle, the pendulum is never at rest, and the small 
quantum fluctuations of the ground state become parametri- 
cally transformed into real, measurable oscillations of the su- 
perconducting phase. According to the Josephson relations, 
these oscillations of the phase produce an oscillatory voltage 
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between the signal line and the ground of our transmission 
line, which propagates as real photons in the transmission 
line. Hence, the Josephson metamaterial acts as an antenna 
broadcasting information about the local vacuum. The math- 
ematical analysis of the DCE in a Josephson metamaterial is 
presented in the Supplementary material. The flux modula- 
tion leads to coupling of modes a(^f + v) and — v)] 
(see Eqs. (S26) and (S49)). 

Our experimental results on the dynamic Casimir photon 
generation in Fig. [2^ display a strong increase in noise power 
around cu TCB when Pd at uid/2-n is enhanced. The increase 
in the noise level is abruptly saturated when a parametric 
oscillation sets in at Pd = —75 dBm. Below the threshold, the 
emitted power spectrum follows closely the predicted form for 
a lossless cavity 
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where l/\M{v)\ 2 is a factor producing parametric gain (in 
S26), x ( w ) i s the electrical susceptibility of the resonator with 
n as the cavity decay rate (S37), and a denotes the magnitude 
of the flux pumping (S14). The theoretical model can be 
extended/generalized to account for the losses in the cavity 
(see Eqs. (S50-S53)). 

The Casimir photons are generated as correlated pairs with 
frequencies uj± — w < j/2±f, with energy conservation satisfied 
uj+ + cj_ — ujd- To study experimentally the correlations 
between the sidebands to± , we use two vector signal ana- 
lyzers triggered simultaneously to collect in-phase and out- 
of-phase time-domain IQ-waveforms at uj± and with homo- 
dyning phases 8±. We calculate the correlations between the 
IQ data and eliminate the amplifier noise by a procedure de- 
scribed in detail in the Supplementary material. The resulting 
correlations, which we denote by (/+/_}, (I+Q-), etc., turn 
out to be directly related to the correlations of the quadra- 
tures x± = (l/2)(a± + a±) and y± = (l/2i)(a± — a±) of 
the outgoing field {a±,aj_} in the two sidebands. For conve- 
nience, the homodyning phases 6 + and of the two analyz- 
ers and the phase of the pump field have been rotated so that 
(x±y^) = (y±x^) = 0. We have verified that all the mea- 
sured quadrature correlations behave according to the theory 
under this rotation. Since this procedure is performed for ev- 
ery point at which we measure correlations, there is no need 
write down the phases 8± any more. To put the correla- 
tions in a compact form, we introduce the covariance matrix 
a(A, B) — \({A, B}}, where {, } denotes the anticommutator 
and A and B are any of the operators x±, y±. The experimen- 
tal results for correlations are displayed in Fig. [2] The time 
dependence of (/+/_) and (Q+Q-) is displayed in Fig. [2]b, 
and our results for the full correlation matrix are illustrated 
in Fig. [§. 

In order to compare the theoretical predictions for Casimir 
photon generation with our results, Fig. 2d displays the 
pumping power dependence of the correlators a(x+,X-), 
a(x+,x+), and their ratio a(x+,x-)/a(x+,x+). When com- 
pared with theory, a(x+,x+) follows closely 0.5 • DCE + | 
(dotted curve) and the covariance cr(x+,x-) agrees with the 
calculated two-mode correlations (dashed curve). The scaled 
squeezing correlator cr(x+, X-)/a(x+,x+) displays the correct 
linear dependence on the pumping drive amplitude Vd but its 
magnitude falls a bit below the theoretical expectation. 

Finally, we have measured what happens when the res- 
onator is detuned from uud/2 while keeping Pd = const. As 
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FIG. 3: a) Emitted noise as a function of v = u — 0Jd/2 
while tuning A = uj — uid/2 measured around Wd/2=5.4 GHz 
at 50 mK. b) Theoretical prediction form Eqs. (S50-S53) 
evaluated using parameters a — 2n x 65 MHz, ke = 2n x 
15 MHz (at 50 mK), ki = 2?r x 9.2 MHz (at 60 mK); in both 
(a) and (b), the color scale for the noise power is cut at 0.3 dB 
in the vicinity of U)d/2. c) Theoretical prediction of classical 
vacuum calculated at a single quantum of classical noise in 
the internal modes, d) Slices with constant detuning A taken 
from (b) are compared with the theoretical predictions from 
Eqs. (S50-S53); parameter values as above. The inset depicts 
the full, symmetric noise pattern measured at smaller pump 
power over the range —100 < A < 100 MHz and —100 < v < 
100 MHz. 



seen in Fig. 3a, we observe a splitting of the noise spectrum 
into two peaks, which are in good agreement with our theo- 
retical results depicted in Fig. 3b at Te = 50 mK and Ti = 60 
mK for the external (ke) and internal (ki) dissipation tem- 
perature, respectively. A classical description of the measured 
spectra as resulting from uncontrollable additional sources of 
noise would utterly fail to explain the structure presented in 
Fig. 3a. As an example, in Fig. 3c we show the effect of a 
single quantum of classical noise in the internal modes: the 
spectrum becomes strongly asymmetric [31] , due to the asym- 
metric dependence of the emitted cavity power on the internal 
temperature in the "lower arm" in Figs. 3a and 3b (see Eqs. 
S50-53). Thus, the symmetric appearance of the pattern in 
Fig. 3a can be regarded as an additional signature of the 
quantum nature of the photon generation in our system. 

The experimental results at three different values of detun- 
ing are plotted and compared with the theoretical predictions 
in Fig. 3d. Based on these fits, or more precisely from the 
symmetry of the lower and upper peaks, we can determine 
that the metamaterial temperature is not much different from 
the measured T = 50 mK. The observed bimodal structure of 
the spectrum with nearly equally strong peaks at v ~ A can 
be considered as direct evidence of the existence of the dy- 
namical Casimir effect. Furthermore, these observations were 
made in a regime where the parametric gain of the driven 
cavity is nearly unity, in clear distinction to regular para- 
metric amplifiers and optical parametric oscillators produc- 
ing squeezed states. Our results pave the way for EPR-type 
of experiments with correlated microwave photon pairs and 
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for quantum information processing with continuous variables 
[291 ED]. 
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